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Abstract 

In the present study, an investigation is carried out to determine 
the effect of soil-rock and rock-rock foundation systems on 
dynamic response of block foundations under vertical mode of 
vibration. The half-space theory is used for the analysis of 
foundation resting on homogeneous soil and rocks. The finite 
element program having transmitting boundaries is considered 
for layered system considering soil-rock and rock-rock 
combinations. The analysis is carried out in details for soil- 
rock and weathered rock-rock systems and the different 
equations are presented for above combinations. The effect of 
top layer thicknesses, shear wave velocity and eccentric 
moments are also simulated. The rock-rock systems considered 
are sandstone, shale and limestone underlain by basalt rock. It 
is interpreted that as the shear wave velocity ratio increase the 
natural frequency increases and the peak displacement 
amplitude decreases. 

Keywords:- Block foundation, Vertical vibration, Soil-rock, 
Rock-rock, 

I. Introduction 

The geology of earth surface is very vast therefore, it is very 
rare to get homogeneous soil in natural state and it can exist in a 
state with a hard rock at shallow depth, consisting of different 
strata with different properties. The two layered systems which 
are commonly available in nature are: (i) soil-rock system and 
(ii) rock-rock system. The dynamic behavior of machine 
foundation placed on non-homogeneous rocks is very complex 
due to bedding foliation plane with varying strength, fissures, 
joints, and faults. The dynamic force on machine foundation is 
mainly due to operation of rotating type machinery. Due to 
complexity of the problem in machine induced vibrations on 
soil-rock and rock-rock system there is a need of research 
addressing these issues. 

The basic mathematical model used for the dynamic 
analysis of machine foundation resting on homogeneous 
medium is a mass-spring-dashpot model. This model is used by 
many researchers [1, 2] to study the dynamic response of 
footing. Many researchers have used an analytical solution 
which considers the foundation soil system as an oscillating 
body resting on semi-infinite, homogeneous, isotropic, elastic 
half-space. Reissner [3] developed an analytical solution by 
using elastic half-space mathematical model for periodic 
vertical displacement at the center of the circular loaded area. 
Sung [4] considered the effects of the contact pressure due to 
the vertical vibration. Lysmer and Kuhlemyer [5] and Veletsos 



and Wei [6] introduced finite element and boundary element 
solutions. 

Many researchers found that natural frequency of the 
foundation soil system depends on shape and size of the 
foundation, depth of embedment, dynamic soil properties, 
nonhomogeneities in the soil, frequency of vibration etc. The 
effect of footing resting on layered soil was studied by some 
researchers [7, 8, 9]. Baidya [10] developed a solution to 
determine the stiffness of the i-th layer of a multilayered 
system and the applicability of this theory have been verified 
by Baidya and Rathi [11] and Baidya et al. [12] from model 
block vibration test results. However, no studies are available 
in the literature to determine the stiffness and damping of 
foundation resting on various soil-rock and rock-rock system. 

In the present paper, an attempt is made to investigate 
the frequency amplitude response of block foundation in 
layered soil-rock and rock -rock system. The size of machine 
foundation block considered for the present study is 

0. 75x0.75x1.0 m (as per IS 5249). The non-homogeneity 
incorporated in the present study is due to the horizontal 
bedding. In the analysis half-space theory is used for 
homogeneous medium (rock or soil). The system is assumed to 
be composite medium for layered combination. The results 
obtained from the analysis of homogeneous medium and rock- 
rock system shows the variation of normalized stiffness, 
damping and amplitude with dimensionless frequency. The 
rock-rock combinations considered in present study are 
sandstone, shale and lime stone underlain by basalt. The 
graphs are plotted to show the variation of dimensionless 
natural frequency and resonant amplitude with shear wave 
velocity ratio for soil-rock and weathered rock-rock systems. 

II. Theoretical Study 

1. Homogeneous medium 

The frequency dependent stiffness and damping of foundation 
resting on the homogeneous strata as halfspace is determined 
by the theory proposed by Veletsos and Verbic [13]. An 
approximate solution is presented for the steady state response 
of rigid disk with mass which is supported at the surface of a 
viscoelastic halfspace and is excited by a harmonically varying 
vertical force (= P,e ,a ", where t = time; co = circular frequency 
of the excitation; and i = %/-!). The response of the system is 
interpreted in terms of the response of an equivalent spring - 
dashpot oscillator with mass. The dynamic characteristic for 
vertical vibration system can be specified by its circular natural 
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frequency (p z ) and by the percent of critical damping (£). These 
quantities can be computed as 



z = < 

2^ (» + <)] 



(2) 



Approximate equation for function l h' can be obtained as 

{H J { RJl-E/H 
The empirical expression given by Kausel and Ushijima [14] 
for the static vertical stiffness using Eqs. 5, 6, 7 and 8 as 
follows: 



where, B, = dimensionless inertia ratio; g*= dimensionless 
parameter; (/?,) = circular natural frequency of the foundation 
considering the halfspace as massless; c_*= damping coefficient 
of the equivalent oscillator; K, = static stiffness of the disk- 
halfspace system; m = mass of foundation; and M'= mass of 
the equivalent oscillator. 

The displacement amplitude of vertically vibrating 
system can be computed as 

(»,). 



1 



+ 4c:\ 



(3) 



Where, (u z ) a = the static displacement produced by the force 

amplitude P z . 

2. Layered medium 

For the analysis of layered medium, top layer is treated as 
stratum and bottom layer is considered as halfspace. The 
dynamic response of layered system is analyzed by using finite 
element program with transmitting boundaries proposed by 
Kausel and Ushijima [14]. 

The FEM program used for analysis is based on the 
formulation and computer program developed by Kausel [15]. 
A three-dimensional axisymmetric finite element model with 
transmitting boundaries was used to model a rigid, embedded 
circular foundation welded to a homogeneous soil deposit of 
finite depth resting upon a much stiffer rock-like material. 
The vertical static stiffness can be derived as 
k 4GR (4) 
" ~ 1-v 

Static stiffness can be expressed in terms of embedment depth 
and strata depth as 

where, f, g, h = Functions such that f(0) = g(0) = h(0) = 1 , 

R/H = Depth ratio, E/R = embedment ratio. 

The expression for stiffness factor '/ can be derived as 

'£)-'♦"■! 

Function for embedment effect can be determined as 



K v = f^{l + 1.28^{l + 0.47^1 + {o. 8 5-0.28^^"j (9) 

where, R = radius of footing, H = depth of soil layer below 

footing, E = embedment depth, 

G = shear modulus, and v = Poisson ratio. 

The equation for dynamic stiffness coefficient can be derived 



k = k v (k'+iac')(l + 2if3) (10) 
where, k v = static stiffness, ft = material damping, a v = 
dimensionless amplitude and k',c' = dimensionless function 
referred to as stiffness coefficients. 

The computer program DYNA 5 which is formulated 
by Novak et al. [16] is used for the present study. This program 
is used to present the dynamic behaviour of block foundation 
as frequency response curves for vertical displacement, 
stiffness, and damping constants for both homogeneous and 
layered system. 

III. Results and Discussions 
1. Homogeneous medium 

The analysis is carried out for homogeneous soil and rocks by 
using the theory proposed by Veletsos and Verbic [13], for 
eccentric moments (m e e, where m e = mass of eccentric rotating 
part in oscillator and e = eccentricity of rotating part of 
oscillator) of 0.278, 0.366 and 0.450 N-m. Four different rocks 
(sandstone, shale, limestone, basalt) are considered for 
analysis. The results so obtained are then plotted in 
dimensionless form. 
Dimensionless amplitude 

X ° = ~tj (11) 
Dimensionless frequency 

_ °> R o 

a - d2) 

where, u z = amplitude of system; v s = shear wave velocity; 
co = angular operating frequency; R„ = radius of foundation 

To normalize the stiffness and damping constant, the 
stiffness and damping values are divided by stiffness of soil at 
zero dimensionless frequency (a„ = 0) 
i _ k 

(13) 



k Soil (a o =0) 
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Where k„ = Normalized stiffness constant, c a = Normalized 
damping constant, and k Soi i(a = 0) = Stiffness of soil at a a = 

The properties of soil and rocks [17] considered for the 
analysis are shown in Table. 1. Figure 1 shows the variation of 
normalized stiffness with dimensionless frequency for soil and 
different rocks respectively. It has been observed that the basalt 
(igneous rock) exhibits high stiffness value as compared to the 
other sedimentary rocks (sandstone, shale and limestone) and 
soil. This is because of the high values of density and the shear 
wave velocity of the basalt as compared to other sedimentary 
rock. 




tin. I \ , .11:. Hi Ml ol normalized sIllIlK"--. '.' Illl illlllellMoliless frequency 




Fig. 2 Variation of normalized damping x\ itli dimensionless tiequene\ 
Table 1 . Properties of soil and rocks 



Type of soil and 


Shear wave velocity 


Unit weight 


rock 


(m/s) 


(N/m 3 ) 


Soil 


185 


16000 


Weathered rock 


1680 


20810 


Sandstone 


1110 


22000 


Shale 


2220 


23000 


Limestone 


2960 


25000 


Basalt 


3700 


26000 



Figure 2 shows the variation of normalized damping 
with dimensionless frequency for rocks and soil respectively. 
As the frequency approaches to zero the damping shows an 
increasing trend. This is because of the conversion of 
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frequency-independent material damping (/?) and frequency 
(co) to the equivalent viscous damping coefficient (c) as 
c = 2j3/co. 

The variation of dimensionless amplitude with 
dimensionless frequency is shown in Figure 3. It is observed 
that the dimensionless resonant amplitude of soil is higher as 
compared to rocks. The dimensionless resonant frequency for 
soils is less than the rocks. Figure 3 can be used to find the 
frequency amplitude response of block foundation of any mass 
and size for homogeneous soils and rocks under vertical 
vibration. 




Fig. 3 Variation of dimensionless amplitude with dimensionless frequency 

2. Layered medium 

The ratio of depth of top layer (H) with width of footing (B) is 
taken as (H/B) = 1.0, 1.5, and 2.0 to investigate the dependency 
of frequency amplitude response on depth of top layer. The 
three different eccentric moment considered in analysis are m P e 
= 0.028, 0.037, and 0.045 kg-m. The results are plotted in the 
form of dimensionless parameters i.e dimensionless frequency 
a , dimensionless amplitude x a etc. The normalized 
stiffness(k ) and damping(c ) of layered system are obtained 
by the following equations 

k - k > 

k BasaIt {a o =0) (15) 

c i 

C " = k Basa „(a,=0) (16) 

Where, k„ and c are normalized parameters for 
stiffness and damping of layered system, k t and c, are 
frequency dependent stiffness and damping constant of layered 
system, k Basa i,(a = 0) is the stiffness of basalt rock at the zero 
dimensionless frequency. 

Soil-rock system and weathered rack-rock system 
The analyses are performed by considering three different 
shear wave velocity ratios (0.8, 0.6 and 0.3). The dynamic 
response of footing are obtained for the above mentioned shear 
wave velocity ratio and H/B ratio. The trend lines are then 
plotted for the graphs between (i) maximum amplitude and 
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shear wave velocity ratio and (ii) natural frequency and shear 
wave velocity ratio. The graphs are extended for the lower 
values of shear wave velocity ratios ( V tI /V s2 <0.5 , where V sl = 
shear wave velocity of top layer, and V s2 = shear wave velocity 
of half-space) by using the trend line equations, to simulate the 
different soil-rock combination and weathered rock -rock 
combination systems. 




Fig.4 Variation of dimensionless natural frequency with shear wave velocity 




Fig. 5 Variation of dimensionless resonant amplitude with shear wave velocity 

Figure 4 shows the variation of dimensionless natural 
frequency with shear wave velocity ratio. The natural frequency 
is found to be constant for all vertical force of excitation for the 
given H/B ratio as nonlinearity is not considered in the analyses. 
It can be seen from the figure that all three trend lines are 
approaching towards the one constant value, which can be 
obtained at the shear wave velocity ratio of 1. Figure 5 shows 
the variation of dimensionless resonant amplitude with shear 
wave velocity for different H/B ratio. The variations of 
dimensionless resonant amplitude with shear wave velocity for 
different H/B ratio are shown by trend line equations. It is found 
that with increase in depth of top layer (H) for the particular 
shear wave velocity ratio the dimensionless resonant amplitude 
value decreases. 



Rock-rock system 

The present study is focused only on non-homogeneity due to 
horizontal bedding plane. Therefore the analysis is carried out 
on the combinations of sedimentary rock underlain by basalt 
(igneous rock) which is most common in geological aspect. 
Sedimentary rocks considered in this study are the sandstone, 
shale, and limestone. 

The variation of normalized stiffness and damping 
with dimensionless frequency (a ) for different H/B ratios and 
for different rock-rock combinations are shown in Figures 6 
and 7 (sandstone underlain by basalt rock, shale underlain by 
basalt rock and limestone underlain by basalt rock) 
respectively. 




Fig.6 Variation of normalized stiffness with dimensionless frequency (rock- 
rock) 




Fig.7 Variation of normalized damping with dimensionless frequcncx (rock- 
rock) 

It is observed that the normalized stiffness value 
decreases with dimensionless frequency (a„). It can be seen 
that with increase in the depth of top layer the normalized 
stiffness decreases i.e sedimentary rock layer up to 
dimensionless frequency value equals to one. It is found from 
the figures that the stiffness value is less for sedimentary rock- 
basalt system than the homogeneous basalt rock system. As the 
frequency approaches to zero the damping shows an increasing 
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trend. The variations of normalized damping with H/B ratios are 
found to be significant for sandstone-basalt combination as 
compared to the other two systems. 

The variations of dimensionless amplitude (x ) with 
dimensionless frequency (a ) are shown in Figures 8 and 9 for 
different sedimentary-basalt rock combinations. It is observed 
that with increase in the depth of top layer, i.e sedimentary rock 
layer, both dimensionless resonant amplitude and dimensionless 
natural frequency decreases. It can be seen that the value of 
dimensionless amplitude is maximum for sandstone-basalt 
system as compared to other two sedimentary rock-basalt 
systems. The variation of maximum dimensionless amplitude 
with H/B ratios is found to be significant for sandstone-basalt 
system as compared to the other two systems. 




Fig. 8 Variation of dimensionless amplitude w ilh dimensionless frequency for 
rock-rock system (Sandstone-Basalt) 




Fig. 9 Variation of dimensionless amplitude w ith dimensionless frequency for 
different rock-rock system 

IV. Conclusion 

In case of block foundation resting on homogeneous systems, it 
is found that both the normalized damping and normalized 
stiffness show higher values in case of homogeneous rocks as 
compared to homogeneous soil. The dimensionless resonant 
amplitude values are found higher and the dimensionless natural 
frequency values are found lower for homogeneous soil as 
compared to homogeneous rock system. 



A detailed procedure has been presented for the 
prediction of dynamic response of block foundation on soil- 
rock combination and weathered rock-rock combination. The 
equations proposed for different H/B ratios in terms of shear 
wave velocity ratio can be used to calculate the dimensionless 
natural frequency and resonant amplitude for soil-rock and 
weathered rock-rock system. It is concluded that with increase 
in H/B ratio, the dimensionless natural frequency and resonant 
amplitude decreases. It is also found that the dimensionless 
resonant amplitude decreases as shear wave velocity ratio 
increases. In the case of rock-rock combination, it is found that 
for all H/B ratios, the resonant amplitude increases and 
dimensionless natural frequency decreases with decrease in the 
shear wave velocity of the top layer. It is observed that the 
normalized stiffness value decreases with increase in both 
dimensionless frequency and H/B ratio. The variation of 
normalized damping with H/B ratios is found significant for 
sandstone-basalt system due to the low values of shear wave 
velocity of top layer. 
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